Observation and Analysis of Resonant Coupling Between Nearly Degenerate Levels of 

the 2^2+ and l^U, States of Ultracold ^^Rbg 

g 9 ^ 

R. Carollo, M. A. Bellos, D. Rahmlow, J. Banerjee, E.E. Eyler, P.L. Gould, and W.C. Stwalley 
Department of Physics, University of Connecticut, Storrs, CT 06269 
(Dated: February 27, 2013) 

We report on the anomalously high line strength of a single rotational level in the ultracold 
photoassociation of two *^Rb atoms to form **Rb2. The v' = 111, J' = 5 level belongs to the 
Hund's case (c) 2(0^) state, which correlates to the Hund's case (a) 2^E^ state. Its strength is 
caused by coupling with a very near-resonant long-range state. The long-range component is the 
energetically degenerate v' = 155, J' = 5 level of the case (c) 2 (Ig) state, correlating to the case (a) 
1 ^Ilg state. The line strength is enhanced by an order of magnitude through this coupling, relative 
to nearby vibrational levels and even to nearby rotational levels of the same vibrational level. This 
enhancement is in addition to the enhancement seen in all J' = 3 and 5 levels of the 2 (0^) state 
due to an / = 4 shape resonance in the a ^Ej state continuum, which alters the distribution of levels 
formed by photoassociation. 

PACS numbers: 33.15.Mt, 33.15.Pw, 33.20.Wr 



I. INTRODUCTION 



Ultracold molecules are currently a topic of much in- 
terest in the atomic and molecular physics community. 
When created through photoassociation (PA) of ultra- 
cold atoms, which is one of the simplest experimen- 
tal techniques, these molecules can be used for high- 
resolution spectroscopy For many other applications, 
there is strong interest in controlling the final state of the 
molecule. In the case of alkali-metal dimers, levels of both 
the X and a'^S+ states are stable enough for study. 
In either state, the vibrational and rotational levels pop- 
ulated are the primary variable that experimenters want 
to control. Several groups have had success in using PA 
to produce molecules in v" = of the X^S^^ state in 

K2JI, Cs2 i, LiCs [3, NaCs 0, KRb 0, and RbCs 
0j Si- own group has previously formed molecules 

in the lowest triplet a'^Sj state of ^^Rb2 in vibrational 
levels v" 32 — 35 ,91 and v" = [l^l , and formation of 
v" = molecules was also previously reported in [ll| . 

There is a significant body of literature regarding en- 
hancement of PA via resonant coupling, and its use is 
becoming important in many experiments. Resonant 
coupling is experimentally useful because the long-range 
state enhances the PA rate, while the short-range state 
enhances decay to desirable (typically deeply-bound) lev- 
els. Our group has previously studied several such reso- 
nant couplings in ^^Rh2 [11 and ^SR^^Rb 0. We have 
also discussed potential applications of these couplings 
and other predicted couplings in the creation of ground- 
state molecules Il3l . Other work using resonant coupling 
includes NaCs |l4l| and Cs2 Here, we demonstrate a 
pathway to form molecules in v" = 18 — 24 of the a^Ti^ 
state via a near-degenerate resonant coupling between 
levels of the Hund's case (c) 2 (0+) state at short range 
and the 2 (1^) state at long range. The 2 (0^) state cor- 
relates with the Hund's case (a) 2 state, while the 
2 (Ig) correlates to the Hund's case (a) 1 ^Hg state. 




FIG. 1. (Color online) Potential energy curves for the 2 (Ojj") 
and 2 (Ig) states as calculated by [3. The levels of the 2 (0^) 
state observed in [l^ are shown, as is the v' = 155 level 
of the 2(lg) state that is energetically near-degenerate with 
v' = 111. 



II. EXPERIMENT 

Our experimental apparatus has been described in de- 
tail in Ref. We begin by loading a magneto-optical 
trap (MOT) with 8 x 10^ '^^Rb atoms at a peak den- 
sity of 10^^ cm~^ and a temperature of 120 fiK. We 
then excite free-to-bound transitions to the 2 (0^) state 
converging to the 5si/2 + 5pi/2 asymptote. It rapidly 
decays radiatively to form metastable molecules in the 
a^^u state, which arc then ionized via 1 + 1 resonance- 
enhanced multiphoton ionization (REMPI) through the 
2 ^T,^ or 3 ^SjJ" states and detected by a discrete dynode 
multiplier. Spectroscopy of these states was previously 
described in Refs. 0, [la]. Molecular ions arc discrimi- 
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natcd from atomic ions and scattered light by time-of- 
flight mass spectrometry. 

As we showed in recent work [13], the 2(0^) state 
supports quasibound vibrational levels behind a barrier 
above the 5s + 5pi/2 limit. As also described in that 
work, the rotational distribution of these levels is affected 
by a ground-state I = 4 shape resonance that enhances 
the photoassociation rate to J' = 3 and 5 beyond the 
strength of lower rotational levels, and well above the 
expected strength for a thermal distribution. All of the 
levels that we have observed are below the 5s + 5p3/2 
asymptote, where there are also many levels belonging 
to other electronic states that correlate to that asymp- 
tote. One such level, v' — 155 (±l|3 of the 2(1^) state 
(shown in Fig. [!}, is nearly energetically degenerate with 
v' = 111 of the 2(0+) state. 

When molecules in the 2 (0+) state decay, they form 



metastable a T,Z molecules in levels 



18 - 24. A 



section of a REMPI spectrum produced from the decay 
of u' = 111 and exhibiting the a"^S+ vibrational level 
spacing is shown in Fig. [2] A similar spectrum showing 
molecules produced by the decay of v' = 107 is shown 
for comparison. While also strong, it has a noticeably re- 
duced signal-to-noise ratio, indicating the v' = 111 level's 
usefulness in spectroscopic applications. Also shown is 
a REMPI spectrum obtained by photoassociating to a 
small satellite peak of the v' = 111, J' = 5 level. This is 
discussed in more detail in Section HVl where the satellite 
peak is marked in the PA spectrum in Fig. |31[b). For a 
pure = state such as the 2 (0+) state, there should 
be no significant hyperfine structure. Nonetheless, de- 
cay products of the satellite peak show nearly the same 



^1]+ state vibrational level distribution as the 



111, 



J' = 5 level main peak, showing that they are closely re- 
lated. This unexpected substructure around v' = 111, 
J' = 5 is one of several indications that the long-range 
2 (Ig) state is coupled with the 2 (0+) state. 



III. MOLECULE PRODUCTION 

A typical PA transition to a molecular level red- 
detuned from the 5s -I- 5pi/2 or 5s -I- 5p3/2 asymptote 
can create trap loss in the MOT of several percent, and 
extremely strong PA lines can exceed 50% loss. By com- 
parison, the blue-detuned PA reported in our previous 
work [loj and in the current experiment [l7| has never 
produced observable trap loss signals. "Blue-detuned" 
PA is used to denote PA to levels that are energetically 
above the atomic asymptote to which they correlate, and 
thus are quasi-bound. Other investigations of this spec- 
tral region by trap loss have also seen no evidence of 



^ The vibrational numbering is somewhat uncertain, as it is based 
on ab initio potentials and no complete experimental assignment 
is known. We believe, however, that the assignment is good to 
within itl vibrational quantum number. 
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FIG. 2. (Color online) Three REMPI spectra of the same 
region, using different PA levels to form the molecules. At 
top (black): a spectrum taken using PA to the resonantly 
coupled v' = 111, J' = 5 level of the 2 (0^) state. At center 
(red): a spectrum using a typical vibrational level, v' = 107, 
J' — 3, which is not resonantly coupled. At bottom (blue): 
the spectrum taken using a small satellite feature of the v' = 
111, J' = 5 line, which is marked in Fig. \MP)- This spectrum 
is the average of three scans and is enlarged by a factor of 10 
for visibility. We believe that this satellite feature represents 
a level of coupled 2 (0^) and 2 (Ij,) character. Vibrational 
quantum numbers of the initial a ''S+ and intermediate 2 '^E+ 



states are assigned as indicated. 



blue-detuned PA [T^. Nevertheless, the strongest line 
in this work, v' = 111, J' = 5, yields REMPI signals of 
275 ions per REMPI shot or more under favorable condi- 
tions. The peak of this transition is clipped and therefore 
this value is actually a lower bound on the actual rate. 
Modeling the clipped line as a pure Lorentzian gives an 
estimated peak ion production rate of 540 ions per 
shot. An undipped spectrum has been scaled to match 
this peak value in Fig. ^h), where the line strength of 
the non-coupled lines (J' = 0-3) is seen to be compara- 
ble to the v' = 110 spectrum in Fig. [Slja). We believe 
this similarity should exist because the ratios of the lines 
that arc not resonantly coupled arc similar, indicating 
that PA is not strongly affected by other factors such as 
the I = 4 ground-state shape resonance, and that the 
resonant coupling affects only J' = 5. 

Using the simple method of Ref . [l3| , we can estimate 
the PA rate leading to this REMPI signal. The number 
of ions measured per REMPI pulse is: 



^Rb^ — A^a-^ionization^d 



(1) 



where Cd is the detector efficiency, Pionization is the prob- 
ability of ionization by the REMPI pulse, and Na is the 
number of molecules in the relevant vibrational level of 
the a^'St state. For a conservative estimate of the PA 
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FIG. 3. (Color online) (a) The rotational spectrum of v' — 
110 of the 2 (0^) state as seen in ion detection from the a ''eJ, 
v' — 20 level, a fairly typical example of a strong PA tran- 
sition to the 2(0^) state, (b) The spectrum of v' = 111 
detected from the a^Ej, v' = 22 level, showing the strongly- 
enhanced J' = 5 line. This spectrum has been scaled based 
on Lorentzian fitting of other spectra, which exhibited clip- 
ping. The resulting J' = 3 line strength is comparable to that 
of v' = 110 in panel (a). At the top of panel (b) is the MOT 
fluorescence trap loss spectrum of the long-range 2 (Ig) state. 
High resolution REMPI and trap loss scans of these coupled 
states are shown in Figs. Ob) and (a), respectively. In both 
panels, the lines marked with an asterisk (*) are hyperfine 
"ghosts" from the lower F — 2 ground-state hyperfine level 

of *5Rb2. 



rate, we will assume a detector efficiency of Cd 
though it may be somewhat less. 

The ionization probability is given by: 



1, al- 



ionization 



1 



-Wt 



(2) 



Here W 



^ is the transition rate determined by the 
photoionization cross section a and the photon flux F per 



unit time. In turn, the flux is F = EX/ {hcirr'^), where 
E is the pulse energy, A is the pulsed laser wavelength, 
and r is the radius of the pulsed beam in the interaction 
region. The pulse is assumed to have constant inten- 
sity, as the laser beam profile is highly non-Gaussian. 
There is little published data on ionization cross sections 
in Rb2, especially for two-photon processes. However, the 
cross section of the upper REMPI photoionization step is 
rate limiting, as the bound-to-bound initial step is likely 
saturated, so we will assume that ionization acts as a 
one-photon process. Using the data from Ref. [l^l and 
allowing for significant deviations due to the different ion- 
ization conditions, we will take a = I^q 5 x 10~^®cm~^. 
This results in Pionization — 0.235, with lower and upper 
bounds of 0.125 and 0.799. The remaining terms are ex- 
perimental parameters, with A ~ 655 nm, E = 5 mJ, and 
r = 1.4 mm. 

With these inputs, and using the modeled peak of 



540, the population is Na = 2300 molecules 



in the detected vibrational level. 



= 22 of the a3E^ 



state, with a range of ^ 675-4300 molecules, 
rate can then be determined using 



Nait) = 



RpAPpcpt 



The PA 



(3) 



Here RpA is the PA rate per atom, r is the time molecules 
spend in the REMPI interaction region after formation, 
and PpcF is the Franck-Condon Factor (FCF) that ap- 
proximates the fraction decaying to an individual triplet 
level. In our system, cold molecules spend ~ 5 ms in the 
REMPI region before their velocity and the acceleration 
of gravity carry them out. The FCF for decay to the 
v" = 22 level is 4.36 x lO'^, as calculated bvLEVEL 
8.0 lU using the a^j:+ potential from Ref. [22|. The 
estimated PA rate is thus 1.1 x 10^ molecules per second, 
with a range between 3.1 x 10^ and 2.0 x 10'' s~^. 

An interesting comparison can be made between the 
PA rate calculated above from ion signals and the rate of 
PA at the same laser frequency implied by the observed 
trap loss signal. This trap loss, as seen in Figs. EJb) and 
SJa), is « 2% at the 2 (0+ ),?;' = Ill position. The trap 
loss is w 4% at the largest peaks of the 2 (Ig) state. For a 
MOT loaded in the presence of an extra loss mechanism, 
the relevant rate equation is 



dN 
~dt 



?'load - (7 + '"pa) 



(4) 



where 7 = 1/t, with r specifying the MOT loading time 
without the PA beam, rioad is the loading rate of the 
MOT without PA, and rpA is the loss rate due to PA. In 
a steady state, the atom number Nq = rjoad/ (7 + ''pa)- 
The measured values of r and iVg are t = 2 ± 1 s and 
A^o = 8 X 10^ atoms. Assuming that the PA laser is 
scanned slowly enough that the steady state is always 
maintained, the PA rate per atom is 



rpA = 7 ( ^ - 1 



(5) 
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Using the 2% value, this gives an estimated rate of rpA = 
1.0 X 10~^ s~^ per atom, with lower and upper bounds 
of 6.8 X 10"'^ and 2.0 x 10^^ s^^ per atom, respectively. 
The total PA rate is 4.1 x 10^ molecules per second, with 
bounds of 2.7 x 10^ and 8.2 x 10^ s'K 

It should be noted that some trap loss is due to 
molecules decaying to free atoms, and if some of these 
are recaptured by the MOT we will underestimate the 
true PA rate. The fraction of 2 (0+), v' ~ 111 molecules 
that decay to free atoms is 70%, which could cause the 
PA rate to be higher by a factor of 3. If we use the high 
estimate for Pionization, the resulting 3.1 x 10^ molecules 
per second estimated from the ion counting rate is in only 
slight disagreement with the estimate based on trap loss 
and possible recapture. 

This in turn is consistent with very strong coupling of 
the long-range 2 (Ig), v' = 155 level and the short-range 
2(0^),!'' = 111, J' = 5 level. Since molecules formed 
at short range appear to account for nearly all photoas- 
sociated molecules, the wavefunction must be strongly 
mixed. 

For further comparison, we can calculate the molecule 
formation rate for the 2 (0^), v' = 110, J' = 5 level, which 
does not benefit from resonant coupling, using the same 
methods. In the spectrum of Fig. Elja), the a^'S^,v" — 
20 level is detected. It has a FCF of 6.60 x lO'^, and 
a peak ion signal size of 46 ions per shot. This gives a 
molecule production rate of 6.0 x 10^ per second (with a 
range from 1.8 x 10^ to 1.1 x 10^), an order of magnitude 
less than the 1.1 x 10 molecules per second from the 
coupled level. 



IV. HYPERFINE STRUCTURE AND 
COUPLING 

To investigate the suspected hyperfine structure men- 
tioned in Sec. |lll we undertook a scries of high-resolution 
PA scans through the J' = 5 line. Each successive scan 
was done at a lower PA intensity, to better show fea- 
tures close to the central line. When aligned, as in Fig. 
|3tb), they show significant, consistent structure around 
the central rotational level. As indicated by dotted lines 
showing some of the strongest features, most of the satel- 
lite features correspond to trap loss of the 2 (1^) state 
shown in Fig. IH^a). This structure cannot be directly 
from decay products of the unperturbed 2 (Ig) state, as 
the FCFs for decay of w' = 155 to the a state are non- 
vanishing only for v" = 38 and 39. These two levels have 
never been observed in REMPI in our apparatus, and are 
believed to be photodissociated quickly by the PA laser. 
In addition, the REMPI spectrum of the marked satel- 
lite feature (displayed in Fig. ^ clearly shows the level 
structure of deeply-bound a^E+ state molecules closely 
matching the spectrum of the strong central peak. Thus, 
as in the REMPI spectrum of the J' = 5 satellite fea- 
ture in Fig. [21 this splitting appears to be rotational and 
hyperfine structure induced by coupling to the v' = 155 
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FIG. 4. (Color online) (a) Hyperfine structure of the v' = 155 
level of the 2 (Ig) state, averaged from six trap loss scans. The 
energy zero point is the position of the «' = 111, J' = 5 level 
of the 2(0^) state. Below the data is a model fit using a 
hyperfine Hamiltonian from Ref. [23] and described in the 
text. The model is vertically offset from the experimental 
data, for clarity, (b) A series of scans at progressively lower 
PA intensities through the = 111, J' = 5 level of the 2 (0+) 
state. The scans are vertically offset and smoothed for clarity. 
To the left of the main peak is the satellite feature whose 
REMPI spectrum is shown in Fig. [2] Both panel (a) and 
panel (b) are on the same scale, and referenced to the position 
of the 2(0+), v' = 111, J' = 5 level at 12798.17(3) cm'^ 



vibrational level of the 2 (1^) state. It is also worth not- 
ing that the signal-to-noise ratio when using as little as 
14 mW of PA power is still quite usable for spectroscopy, 
and is produced at a PA laser power far below the 500 
mW to 1000 mW typically used in our experiment. 

As discussed above, our trap loss scans of the v' = 155 
level show significant rotational and hyperfine structure, 
though it is not well-resolved. We believe the linewidth 
is not a result of our scan speed, but is due to some 
combination of natural linewidth and broadening due to 
the high laser power needed to produce observable trap 
loss. An averaged scan, compiled from six high-resolution 
scans, is shown in Fig. [H^a). This spectrum is difficult 
to assign a priori, as the rotational and hyperfine level 
spacings overlap, particularly for low J levels. In order to 
gain insight into this structure, we model these data with 
a simulated spectrum using a Hamiltonian described in 
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Ref. [13, 



= A,{i + fjit + 3)2) + ^ (^F - J- if . (6) 
This can be rewritten as 

n = A,, [i + + 'if) + {F{F + 1) + 2 
+ /(/ + 1) - 2(j)n - F+I^ - F^I+ - 2(t)i + 2ni) (7) 

where Ay is the hyperfine couphng term, arising primar- 
ily from the Fermi contact interaction. By is the 
rotational constant, i, (f>, and fl are the projections of 
I, F, and J on the internuclear axis, respectively, and i] 
is a fitting constant for the term that is quadratic in i 
(this term is added to allow for variation in Ay with i). 
Using the parameters Ay = 2.8 x 10^^ cm^^, rj = 8, and 
By = 0.00095 cm~^, the model qualitatively reproduces 
the observed spectrum, as is shown in Fig. SJ^a). One 
minor issue with this fit is that it requires a rotational 
constant By that is smaller than the By = 0.0012 cm~^ 
calculated from our potentials. The calculated By values 
from other potentials have been relatively accurate, al- 
though they display a slight tendency to overestimate the 
By values compared to experiment [13 • The model 
also does not include weighting of the incoming partial 
waves, which arc affected by the thermal population of 
the MOT as well as the I = 4 shape resonance. Note 
that we have tried to model neither the hyperfine struc- 
ture of the v' = 111 level of the 2(0+) state nor the 
coupling between the two hyperfine-split J' — 5 levels 
(2(0+) ~ 2(1,)). 

From the spectrum in Fig. \M.a), we measure the spac- 
ing between v' = 111, J' = 5 of the 2(0+) state and 
v' = 155, J' = 5 of the 2 (Ig) state to be A « 1.67 x 10"^ 
cm~^. We believe that levels other than J' = 5 of the 
2 (0+) state do not contribute, as the other J levels are 
energetically much further away. If we know the unper- 
turbed spacing of the levels and can measure a shift, we 
can est imate the stre ngth of the coupling interaction via 
A = ^4|i/i2p + Ag, where Aq is the unperturbed spac- 
ing and |-ffi2| is the interaction term of the Hamiltonian. 

By fitting the rotational progression of = 111, J' = 
5, using a simple Ej = ByJ(J+l) model both with and 



without the J' = 5 level, we find a shift of 7 x 10 ^ cm ^. 
This is within the 1.3 x 10"^ cm"^ FWHM linewidth, 
but can still help establish an estimate of the coupling 
strength. Assuming the shift is symmetric with the other 
couphng partner, we find \Hi2\ « 3.4 x 10"^ cm~^, an 
extremely small value. 

Since the two states that are involved in this cou- 
pling have ri = and $7 = 1, the coupling must be 
the result of an inhomogeneous perturbation. There are 
two primary causes of such perturbations — the non-Born- 
Oppenheimer S'-uncoupling and the L-uncoupling oper- 
ators [13] . The S'-uncoupling operator has selection rules 
AS" = 0, Af7 = AE = ±1, and AA = 0. However, 
this operator only couples f2 components of the same 
electronic state multiplet. Additionally, our Hund's case 
(a) 2 ^I]+ and 1 ^Hg states have A = and A = 1, re- 
spectively, ruling out an S'-uncoupling mediated pertur- 
bation. The i-uncoupling operator has selection rules 
An = AA — ±1 and AS = 0. This operator, having 
the form (l/2fiR'^) (J+L" + J"L+), can couple differ- 
ent electronic states. The observed resonant coupling is 
thus likely the result of an L-uncoupling mediated per- 
turbation. 



V. CONCLUSION 

We have shown that there is resonant coupling between 
a pair of J' = 5 levels in the 2^S+ and l^Hg states. 
This coupling causes an order-of-magnitudc increase in 
the production of a "^5]+ state molecules, compared with 
nearby vibrational levels of the 2 ^E+ state, and yields 
an approximate PA rate of 5.4 x 10^ molecules per sec- 
ond. This coupling provides a strong pathway for creat- 
ing deeply bound a'^S+ state molecules. As it connects 
high-u levels of the a'^S+ state (through the potential of 
the long-range component) with more deeply bound lev- 
els (through the potential of the short-range component), 
it can also provide an experimental pathway for molecule 
transfer. 
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